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Abstract. Dynamical coupled-channel approaches are a widely used tool in hadronic physics that allow 
to analyze different reactions and partial waves in a consistent way. In such approaches the basic interac- 
tions are derived within an effective Lagrangian framework and the resulting pseudo-potentials are then 
unitarized in a coupled-channel scattering equation. We propose a scheme that allows for a solution of 
the arising integral equation in discretized momentum space for periodic as well as anti-periodic boundary 
conditions. This permits to study finite size effects as they appear in lattice QCD simulations. The new 
formalism, at this stage with a restriction to 5- waves, is applied to coupled-channel models for the (t(600), 
/o(980), and ao(980) mesons, and also for the yl(1405) baryon. Lattice spectra are predicted. 

PACS. ll.80.Gw Multichannel scattering - 12.38.Gc Lattice QCD calculations - 14.20.Gk Baryon reso- 
nances - 14.40.-n Hadrons, properties of mesons - 24.10.Eq Coupled-channel and distorted- wave models 



1 Introduction 

At low energies chiral perturbation theory has allowed 
for a comprehensive understanding of the strong interac- 
tions. At higher energies, where the chiral expansion starts 
to break down, rich spectra of excited states have been 
found in the meson-meson and meson-baryon sectors ex- 
perimentally, see e.g. Ref. [T]. Several analysis tools have 
been developed to disentangle the partial wave content 
and to pin down the resonance spectrum. Among those 
tools are dynamical coupled-channel models [UlSllllISUSlEl 
[Sll^irO]. characterized by interactions driven by hadron 
exchange, which are derived from effective Lagrangians. 
Such coupled-channel models usually fulfill two-body uni- 
tarity as well as some requirements of three-body unitar- 
ity and crossing symmetry. They are widely used for data 
analyses nowadays, for example by the Jiilich [51I31I31E1E1 
[7], EBAC [S] and DMT [TU] groups in their investigations 
of the TriV system in the second and third resonance re- 
gion. 

In recent years information about the excited hadron 
spectrum is becoming available from the rapidly evolving 
field of lattice gauge theory. The masses of meson and 
baryon ground states could be determined at the few per- 
cent level [TT] and close to physical pion masses. With re- 
gard to excited hadrons, signals of a rich spectrum could 
be found recently [T^[T^lfT41l 1 5 lITB] in calculations with 
pion masses > 300 MeV or in the quenched approxima- 
tion. Of course, those excited states will be able to decay 
on the lattice, once the employed pion mass approaches 



the physical value. Then, resonance signals on the lattice 
manifest themselves in "avoided level crossing"; for nar- 
row resonances, this signal is rather clear, but for broader 
resonances the avoided crossing becomes smeared out [TTj 

In the one-channel- and also two-channel case, and 
given sufficient precision for the levels, the phase shift can 
be extracted unambiguously and without any model input 
by using Liischer's formalism [ropDPT|[^[^Pi| . How- 
ever, the precision of the lattice data is limited. Moreover, 
at higher energies often many more physical channels are 
open as it is the case, e.g., for the ttN interaction in the 
second and third resonance region. Thus, model input may 
be required to extract the physical properties and the res- 
onance spectrum, in particular for broad resonances. 

The physically realistic situation of multiple open de- 
cay channels is most conveniently organized in a coupled- 
channel scheme. Note that this is of particular importance 
in the study of finite size effects as discussed in detail, e.g., 
in Ref. Threshold openings show the same avoided 
level crossing as resonances, and it is, thus, necessary to 
take the physically allowed channels consistently into ac- 
count, at least if they are known to couple strongly to the 
system of interest. 

In Refs. [2^125] a coupled-channel fit to synthetic lat- 
tice data with given errors has been carried out, stabilizing 
the extraction of phase shifts and pole positions by using 
Unitarized Chiral Perturbation Theory to expand the po- 
tential. The feasibility of this procedure has been shown 
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and the error propagation from the lattice data to pole 
positions and phase shifts could be quantified. 

In a similar way, from the point of view of hadron ex- 
change approaches, the energy levels provided by lattice 
calculations can be treated as additional data that enter 
the combined analysis of hadronic reactions. Thus, the mo- 
tivation of this study is to reformulate the formalism of the 
dynamical coupled-channel approaches in a manner that is 
suitable for a prediction of the discrete energy levels on a 
finite-size lattice. In particular, we find that the discretiza- 
tion of the multi-channel Lippman-Schwinger equation in 
the momentum space allows one to achieve the goal with a 
minimum effort. The paper contains a detailed discussion 
of such a discretization procedure. 

The proposed framework can be used to produce the 
volume-dependent energy spectra for different physical sys- 
tems. We are aware that present day lattice simulations 
do not provide more than a few data points for a limited 
number of volumes. Still, we would like to stress again 
that, as shown in Ref. assuming a simple and suffi- 
ciently general parameterization for the chiral potentials, 
it is possible to extract the resonance parameters from a 
fit of the predicted energy levels to reasonably small data 
sets. 

To demonstrate the discretization of dynamical coupled- 
channel approaches, in this study we utilize hadron-ex- 
change potential models that have been developed to de- 
scribe the meson-meson interaction in the strangeness S = 
sector p51[?7] and the meson-baryon interaction for S = 
— 1 as, e.g., in Refs. PSII^ . The meson-baryon sector with 
S = is left to future studies. 

In the aforementioned models, the interaction is given 
by S-, t-, and u-channel exchange of hadrons, which con- 
nect all partial waves of the scattering amplitude. Fur- 
thermore, using SU(3) symmetry, one can also relate dif- 
ferent reactions in the same framework, as recently demon- 
strated [5]. Thus, hadron exchange puts strong constraints 
on the coupled-channel amplitudes and provides at the 
same time a structured background from t— and u-chan- 
nel processes. Together with a minimal set of s— channel 
processes ("genuine" resonances) this allows for a reliable 
extraction of poles and residues from the analytic contin- 
uation of the amplitude 

The considered exchange processes constitute the in- 
teraction potential, V, which is then projected to partial 
waves that represent the basis in which the scattering 
equation is solved. Omitting channel indices, the scatter- 
ing (Lippmann-Schwinger (LS)) equation in the center-of- 
mass frame reads 



OO 

+ / dqq^V{q",q) 



- Eaiq) - Ebiq) + ie 



-T{q,q') , 



(1) 



where q" = \q"\ {q' = \q'\) is the modulus of the outgo- 
ing (incoming) three-momentum that may be on- or off- 



and Eb = \/m^^+q^ are the on-mass shell energies of the 
intermediate particles a and b. Note that Eq. ^ is for- 
mulated in the partial wave basis, i.e. the amplitude only 
depends on the modulus of the incoming, outgoing, and 
intermediate particle momenta. The angular dependence 
of the full T-matrix is, for meson-baryon interaction, pro- 
vided by the Wigner d^^y (??)-functions in the partial wave 
decomposition [28j . where •& is the scattering angle and 
A (A') is the helicity of the incoming (outgoing) baryon. 
For meson-meson interaction, the angular dependence re- 
duces to the standard Legendre polynomials Pj {cos -d). 
The solution T of the LS equation allows to calculate the 
observables. Note that the pseudo-potential V that ap- 
pears on-shell and also half off-shell in Eq. ^ is fixed 
from the underlying interaction Lagrangian. In the follow- 
ing, we abbreviate the denominator in Eq. ([T]) as ^/s— Ei^t , 

i.e., Eint = Ea{q)+Eb{q). 

Dynamical coupled-channel models of the ttN interac- 
tion often take into account intermediate and final ttttN 
states [HISIEIIUIMZIIH] I though usually only in an effective 
way. Still, this complicates the structure of Eq. dH) [5]. 
Thus, for the present study we concentrate on physical 
problems where such three-body channels play only a mi- 
nor role. The two interaction models we consider here [261 
[25] do not contain any effective two-meson states. Also, 
note that the spherical symmetry on the lattice is broken, 
and partial waves mix; however, this effect is small for S- 
waves which mix only with waves of an orbital angular 
momentum L > 4 [5D]. In this study, we concentrate on 
the formal development, while the issue of partial wave 
mixing and multi-meson intermediate states is left for fur- 
ther investigations. 



2 Discretization of the momentum space 

A standard way of solving the integral equation ^ for T 
is to discretize the integral and invert a large matrix in 
the space of off-shell momenta. We will show in the fol- 
lowing that with a few minor changes this scheme can be 
adapted to the discretized momentum space that corre- 
sponds to finite-size lattices as used in lattice gauge cal- 
culations. Explicitly, the integral term in the one-channel 
case is written in terms of the matrices V, G, and T 



^^^,Viqy)nq,q') 
u 

Vij = V{qi,q.j) , Tij = 
G = diag 



V^-Ei' 



^ VGT, 

T{qt,qj) , 



(2) 



where the qiS are the sampling points (i,j = I,-- - ,7i) 
of the integration with the associated integration weights 
Wi, and Ei is i?int evaluated at qi. The multi-channel case 
is discussed at the end of this section. The LS in its dis- 
cretized form, 



shell, y/s is the scattering energy, and Ea ~ \J vri\ + 9^ 



VGT 



(3) 
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reads component-wise for k ~ 1, 



V,; 



fc=l 



V~s-Ek 



-kj ■ 



(4) 



The solution T can now be obtained by matrix inversion, 
i.e. 



(1 - VG)-^V . 



(5) 



For the calculation of phase shifts or observables one needs 
the solution of the reaction amplitude T at the on-shell 
momenta. To obtain the on-shcU result one can augment 
the matrices above by adding the on-shell momentum to 
the sampling points, i.e., Qn+i = ^on, where qon fulfills 

= V?7ia + 9on + V'^b + lln- Then, T„+i,„+i is the 
on-shell amplitude and Ti_„_|_i (i = 1, • • • ,n) are the half 
off-shell amplitudes. Note that Gn+i,n+i = 0. Further- 
more, the two-body singularity occurring at the on-shell 
energy in the LS equation ^ needs to be treated numer- 
ically. One possibility consists in rotating the sampling 
points into the complex plane [30| . Another commonly 
used option is the Haftel-Tabakin scheme [31] in which 
the singular term in the evaluation of the integral in Eq. 
^ is subtracted numerically and added analytically at 
the on-shell point. 

For the discretization on the lattice, none of these two 
schemes needs to be employed, as we will see in the follow- 
ing. The lattice allows only for discrete momenta (periodic 
boundary conditions) which, for a simple cubic lattice, re- 
sults in the substitution 



L3 
27r 

— n, n G 
L 



(6) 



To apply the discretization to dynamical coupled-channel 
models, the LS equation ^ is rewritten as a three-di- 
mensional integral and discrctizcd according to Eq. ©. 
The three-dimensional summation can be further simpli- 
fied by considering the multiplicities of the i-th neigh- 
bors for the simple cubic lattice with periodic bound- 
ary conditions, labeled as (P) in the following. It can be 
shown that these multiplicities are given by the -series 
[^(P)(i = 0,1,2,- • •) = l,6,12,---],seeRcf. [22]. The t?(P) 
series is given by the coefficients of the Taylor expansion 
around a; = of the function g'^^\x) ~ ["i^sCO, xyf ^ 



i=0 



(7) 



where is the elhptic d function [221, 

oo 



(8) 



One obtains the scattering equation in discretized momen- 
tum space. 



T^^\q",q')^V{q\q') 

Viq",q,)T(PHq. 



i=0 



Ea{qi) - Eb{qi)' 



ii -r 



27T 

T 



(9) 



Please note that the sampling points g, of the summation 
in Eq. Q are different from those in Eq. ^ {q.i is the dis- 
tance to the i-th neighbors). Note also that all quantities 
in this equation are real. Furthermore, T*^^^ will now have 
an infinite tower of poles, corresponding to the discrete 
spectrum of a system in a finite volume. 

Finally, we can formulate convenient replacement rules 
to obtain Eq. Q from existing dynamical coupled-channel 
models of the form of Eq. © as given in Refs. [2 1 [5 1 [i l [5 1 [ ^ 
[26.,27..28.,29j. Comparing Eq. (gj) with Eq. (O, we imme- 
diately obtain for the integration measure q^ Wi and for 
the integration sampling points at q.i 



2^2^(P)(i) 



2tt 



L3 

/i, i = 0, 1,- 



(10) 



i.e. the first sampling point is at = with i?'^^)(0) = 
1. Note that if V and T do not depend on q, as it is 
the case, e.g., in the chiral unitary approaches considered 
in [33''341 13!51IHM57[IH51[^ . those quantities factorize from 
the integral in Eq. ^ so that the integral equation ([1]) 
reduces directly to an algebraic equation. The discretized 
amplitude can then can be rewritten, up to order 
effects, in terms of the Liischer zeta function Zqq within a 
/^-matrix approach and for coupled channels as developed 
in Refs. [221I251IM] . The influence of some suppressed 
effects has been recently discussed in Ref. [23]. Also in 
the present formalism, with a non-factorizing scattering 
equation, it would be possible to rewrite the discretized 
LS equation ^ in terms of Zqq (up to effects of order 

Hybrid boundary conditions [101HT1H2] were introduced 
to distinguish scattering states from tightly bound quark- 
antiquark systems. Similarly, as proposed in Refs. [231 
[2i] . twisted boundary conditions provide the possibility 
to change thresholds in lattice gauge calculations. This 
provides a unique opportunity to study the nature of res- 
onances that lie close to a threshold [i51l44U45ll46llTf] like, 
for example, the /o(980) with regard to the KK thresh- 
old, because the twisting moves the threshold while the 
resonance stays put. We realize that it could be quite chaL 
longing to implement this idea (including twisting for the 
sea quarks) in present-day lattice simulations. 

Twisted boundary conditions, e.g. for the strange quark, 
lead to a change in the summation of the lattice momenta 
of the KK channel EHl, 



A;— — OO 



(11) 
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where is the twisting angle that can be chosen to vary 
the KK threshold. For the present study of dynamical 
coupled-channel approaches, it is possible to formulate 
a scattering equation similar to Eq. ([5]), but with maxi- 
mally twisted (i.e. anti-periodic) boundary conditions — 
(tt, tt, tt). For this value of 0, the summation exhibits again 
a high symmetry: compared to the periodic case, the ori- 
gin of the summation is simply shifted by (1/2, 1/2, 1/2), 
i.e. it is in the center of the cubic unit cell. The mul- 
tiplicity of the next neighbors is 8, there are 24 next- 
to-next neighbors and so on. In general, the multiplic- 
ity of the j-th neighbors is given as 8 times the number 
of ordered ways of writing n as the sum of 3 triangular 
numbers. We call this multiplicity d^^^ in the following 
and (i = 0, 1, 2, • • • ) = 8, 24, 24, • • • . It can be shown 
that this series is given by the coefficients of the Tay- 
lor expansion [cf. Eq. ([7])] around x = of the function 
g^^\x) = [i92(0, V^)]^a;"3/^ where ??2 is the elliptic d 
function [321. 



^2(0,VJ) = 2x1/8 



CSO 

^ ^fc(fc+l)/2 
fc=0 



(12) 



With the new multiplicities 'd^^\ and taking into account 
that the sampling points qi of the summation also change, 
the scattering equation with anti-periodic boundary con- 
ditions takes the form 

T<«'(9".9') = V(5",5') 

L-^ V-s - Ea(qi) - Eb[qi) 



with redistributed sampling points qi 



27r 

T 



a + 3 



0,1, 



(14) 



As before, qi is the distance to the i-th neighbors. Like 
for the periodic case, simple substitution rules with re- 
spect to the continuum formalism can be formulated. The 
integration weights change according to 



27r2,?(A)(z) 



L3 



(15) 



and the distribution of the qi is given by Eq. (fH)) for the 
simple cubic lattice with anti-periodic boundary condi- 
tions. 

So far, the discretization has been shown only for the 
one-channel case. The multi-channel case leads to no fur- 
ther complications and the formalism applies to it without 
changes. One just has to take into account that the quan- 
tities appearing in Eq. ^ depend also on the channel, 
for example Vy — >■ V^'', where fj, and v characterize the 
outgoing and incoming channels. The sum runs now over 
both the momenta and the intermediate channels. The 
matrix G, cf. Eq. ([2]), is diagonal in channel space, too. 
If anti-periodic boundary conditions are imposed for, e.g., 
the strange quark in meson-meson scattering P51l24j . the 



substitution rules of Eq. (flSl) are applied to the KK chan- 
nel, while for the tttt channel the periodic rules of Eq. ([TU)) 
are used. Note that in this case one has transition po- 
tentials V'^'^{qi,qj) with periodic {qi, fi = tttt) as well as 
anti-periodic (qj, v = KK) sampling points together. 



3 Results and Conclusions 

In lattice calculations, the energy levels E{L) can be ob- 
tained as a function of the box size L of the simple cubic 
lattice, which allows to apply the Liischer formalism \^ 
to extract phase shifts. For the extraction of poles and 
phase shifts in the multi-channel case, see also Refs. [231 
[23]. Here, we simply predict the levels E{L) applying the 
formalism developed in the previous section. The levels 
E{L) are given by the singularities of the scattering equa- 
tion © or i.e. by the zeros of 



dct(l - VG) = 



(16) 



where V and G are now given by applying the substitution 
rules of Eq. pU]) or Eqs. ([TJ]) and ([T^ to the quantities 
of Eq. ^ , depending on whether periodic or anti-periodic 
boundary conditions arc employed. Note that in a practi- 
cal calculation the sums in Eqs. © and need to be 
cut off at some imax which is chosen in a way that the 
regularization schemes render contributions from i > imax 
negligibly small (this can always be achieved). For the po- 
tential models studied in the present paper imax had to be 
chosen in such a way that w 2.5 - 3 GcV. 

We consider the two dynamical coupled-channel mod- 
els mentioned in the Introduction. Specifically, we utilize 
the interaction potential developed in Refs. [^I?7] for the 
/ = 0, 1, S = meson-meson sector for which the coupling 
between the tttt (or ttt]) and KK channels is taken into ac- 
count. The cr = /o(600), /o(980), and ao(980) resonances 
all appear dynamically generated, i.e. without the need 
of corresponding s-channel exchanges. On the other hand, 
the model does include also a genuine resonance in form 
of an s-channel pole diagram, tentatively called /o(1400) 
in Ref. [55], whose bare mass is at 1520 MeV. 

In the / = 0, S = —1 meson-baryon sector, the dynam- 
ical coupled-channel model of Refs. [^[2U] is considered. 
It includes the channels ttA, ttE, and KN. Like in the 
meson-meson sector, the attraction in coupled channels 



^ Note also that the inverse of qt^^^ corresponds to distances 
of around ~ 0.07 — 0.08 fm. At these distances the form fac- 
tors that provide the regularization render the high-g contri- 
butions so small tliat the results practically do not change. 
Consequently, if the lattice spacing a is taken less or equal to 
the above values (which are common values in present-day lat- 
tice simulations), the finite-a artifacts are washed out and the 
limit a — >■ can be safely performed, as was done in the present 
paper from the beginning. Still, a comprehensive treatment of 
systematic uncertainties tied to the finite spacing would require 
to address the actual lattice action that generates the levels, 
but that is beyond of what can be possibly done in the present 
framework. 
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Fig. 1. Spectrum E{L) for the 7 = 0, = meson-meson 
sector [cr(600), /o(980)] for periodic boundary conditions. Solid 
lines: results for the dynamical coupled-channel approach of 
Refs. [26II27| . Dashed lines: For comparison, the results for the 
chiral unitary approach of Ref. [23], see also Refs. |23ll24j . Dot- 
ted (dash-dotted) lines: free tttv (KK) levels. 



is so strong that the yl(1405) appears dynamically gener- 
ated. In this approach, also the channels KA, K*N, and 
K*A are included. However, these channels contribute 
only effectively to the direct KN interaction in form of 
box diagrams Below the respective thresholds [E < 
1700 MeV) of these intermediate states, the discretizcd 
sum of Eq. ^ converges towards the integral in the large 
L limit and, therefore, at this stage of the exploratory in- 
vestigation we refrain from an explicit discretization of the 
integrals related to these box diagrams. 

The calculated levels, using periodic boundary condi- 
tions, are shown by the solid lines in Figs. [1] [51 and [31 to- 
gether with the respective free levels of the involved chan- 
nels. In case of the meson-meson interaction the results 
exhibit a characteristic avoided level crossing at around 
980 MeV (more pronounced for 7 = than for 7=1). 
In general, avoided level crossing is a signal for a reso- 
nance, but the avoided crossing is washed out for broader 
resonances. This applies for example to the broad /o(600) 
resonance as is obvious from Fig. [H But it is also the case 
for the genuine /o(1400) resonance which acquires a signif- 
icant width due to renormalization in the course of solving 
the Lippmann-Schwinger equation ([T]) . Its pole is finally at 
1346±i 249 MeV, cf. Table V in Ref. [Hj, and thus it is al- 
most as broad as the /o(600). In both cases no direct sign 
of the presence of these resonances is visible in the levels 
[cf. Fig. [Tj. Nevertheless, in case of the /o(600) the pole 
position can be reconstructed from the levels (with large 
errors) by expanding a general pseudo-potential in powers 
of the scattering energy, as recently shown in Ref. [H] . In 
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Fig. 2. Spectrum E(L) for the / = 1, = 0+ meson- 
meson sector [ao(980)] for periodic boundary conditions. Solid 
lines: results for for the dynamical coupled-channel approach 
of Refs. [26II27| . Dashed lines: For comparison, the results for 
the chiral unitary approach of Ref. [33]. Dotted (dash-dotted) 
lines: free Trrj (KK) levels. 
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Fig. 3. Spectrum E{L) of the I = 0, S = -1, = ^ meson- 
baryon sector [yl(1405)] for periodic boundary conditions. Solid 
lines: results for the KN Jiilich model of hadron exchange [281 
I29| . Dashed lines: For comparison, the first two levels as ob- 
tained from the chiral unitary approach of Ref. [33 • Dotted 
(dash-dotted) lines: free irS (KN) levels. 
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the future, we plan to address this (inverse) problem of 
how to reconstruct resonances from a given set of lattice 
data, in the present framework of hadron exchange. 

Avoided level crossing is also present at particle thresh- 
olds. Resonances close to thresholds, like the /o(980) that 
is close to the KK threshold, require, thus, special atten- 
tion. The avoided crossing, visible in Fig. [U is therefore 
not a clear signal for the /o(980), but necessitates ad- 
ditional analysis. The inverse problem of reconstructing 
phase shifts and resonance poles from lattice levels, in such 
circumstances, has been addressed recently in Ref. [23], see 
also Ref. gS]. 

As mentioned above, the channel space in the / = 
sector comprises the tttt and the KK channels. The 47r 
channel has been neglected as there are no phenomeno- 
logical indications by the PDG [SD] that require its inclu- 
sion at the energies of interest. Still, at energies beyond 
the /o, the analysis could be refined by including these 
states, e.g. in terms of quasi-particle channels such as crtr, 
pp, • • • , as carried out, e.g. in Ref. |51j. In any case, the 
level structure of 4-body states is considerably more com- 
plicated than the one of two-body states considered here, 
and beyond the scope of this study. 

In the meson-baryon sector, cf. Fig. [31 the level struc- 
ture below the KN threshold is interesting because of 
a conjectured two-pole structure of the yl(1405) 
Indeed, also the hadron-exchange approach of Ref. pS] , 
considered here, predicts two poles below the KN thresh- 
old [55]. Note also the recent work of Ref. [S2| in which, 
for the first time, the yl(1405) could be isolated in a lattice 
QCD calculation near the physical pion mass. 

For the 7 = 5 = meson-meson case, we also ap- 
ply the formalism for anti-periodic boimdary conditions 
developed in the previous section. Corresponding result 
are presented in Fig. [3J The bending of the levels at the 
KK threshold from avoided level crossing disappears with 
anti-periodic boundary conditions as indeed clearly visi- 
ble in the figure. This demonstrates again the usefulness of 
anti-periodic boundary conditions to study near-threshold 
resonances like the /o(980) discussed above [^1^ . 

For comparison, in Figs. [TJ [51 and [H we also show the 
spectra calculated in Refs. [5H1I5^ , using the chiral unitary 
approach from Ref. [33]. As seen in the figures, the lev- 
els are similar to the result based on the hadron-exchange 
interaction but not identical, see Fig. [T] Indeed, also the 
respective phase shifts in the continuum case show no- 
ticeable differences [2^133] . The further interpretation of 
the calculated spectra and observed small discrepancies 
between different approaches are left for future studies. 

In conclusion, dynamical coupled-channel models can 
be modified so that the calculation of lattice spectra be- 
comes possible, in agreement with the Liischer formalism 
up to effects of order . This framework can be used to 
predict lattice spectra, as done here, or to analyze lattice 
spectra once lattice data will become available. Dynami- 
cal coupled channel approaches, which respect analyticity, 
unitarity and other general requirements of the 5-matrix, 
can thus provide the opportunity to analyze and interpret 
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Fig. 4. Same as Fig. [T] (meson-meson, / = 0) but with anti- 
periodic boundary conditions (6 = tt) in the KK channel, for 
the dynamical coupled-channel model [26II27| (solid lines). For 
comparison, the result with anti-periodic boundary conditions 
for the chiral unitary approach of Ref. [33] is shown (thick 
dash-dotted lines). Also, the result with periodic boundary 
conditions of Fig. [T] is again shown (dashed lines). Dotted (thin 
dash-dotted) lines: free tttt (free anti-periodic KK) levels. 

both experimental data and lattice data within one single 
approach. 
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